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Abstract. We present interferometric observations of 
CO (1-0) emission in the elliptical galaxy NGC 759 with an 
angular resolution of 3 '.'1x2 '.'3 (990 x 735 pc at a distance 
of 66Mpc). NGC 759 contains 2.4 x IO^Mq of molecular 
gas. Most of the gas is confined to a small circumnuclear 
ring with a radius of 650 pc with an inclination of 40° . The 
maximum gas surface density in the ring is 750 M© pc^^. 
Although this value is very high, it is always less than 
or comparable to the critical gas surface density for large 
scale gravitational instabilities. The CO J=2— 1/J=1— 
line ratio is low (0.4), implying sub-thermal excitation. 
This low ratio is consistent with a two-component molec- 
ular gas, consisting of a cold and dense phase, containing 
most of the mass, and a warm and diffuse gas compo- 
nent, dominating the CO emission. Compared to galaxies 
of similar gas surface densities, NGC 759 is underluminous 
in LpiR with respect to its molecular gas mass, suggest- 
ing that the star forming efficiency in NGC 759 is low. We 
discuss the possibility that the molecular gas and current 
star formation activity in NGC 759 could be signatures of 
a late stage of a merging between two gas-rich disk galax- 
ies. We use a mass model of the underlying galaxy which 
is applicable to spherical galaxies with an r^' "^-luminosity 
profile when interpreting our CO data. This leads to more 
modest estimates of the molecular gas mass fraction and 
the surface gas density than would have been derived using 
simpler models, suggesting that many of the spectacular 
molecular gas properties of ultraluminous FIR galaxies, 
which could be described by similar mass distributions, 
may have to be revised. 

Key words: ISM: molecules - Galaxies: elliptical and 
lenticular, cD - Galaxies: individual: NGC759 ~ Galaxies: 
ISM - Galaxies: kinematics and dynamics 
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1. Introduction 

Faber & Gallagher (1976) found the apparent absence of 
an interstellar medium (ISM) in elliptical galaxies surpris- 
ing since mass loss from evolved stars should contribute as 
much as 10^ - 10^" M© of gas to the ISM during a Hubble 
time. Today we know that ellipticals do contain an ISM, 
albeit very different from that found in spiral galaxies. 
The most massive ISM component is a hot X-ray radi- 
ating halo, with a typical mass of 10^ — 10^° Mq (Fab- 
biano 1990). Small amounts of photoionized gas (typically 
lO'^ — 10** Mq) can be found in about 60% of all ellipticals 
(Caldwell 1984, Phillips et al. 1986). In addition, ~ 45% 
of a magnitude limited sample of ellipticals were detected 
at both 60 and 100 ^m by IRAS (Jura et al. 1987, Knapp 
et al. 1989), suggesting the presence of a significant dust 
component. Atomic hydrogen gas is seen in ^ 15% of all 
ellipticals (Knapp et al. 1985). A similar detection rate has 
been found for molecular gas, seen through CO emission 
(Lees et al. 1991, Wiklind et al. 1995). 

Elliptical galaxies have less molecular gas per blue 
luminosity than spirals and a much larger dispersion in 
their L-g/M^i^ ratios. Wiklind et al. (1995) found that the 
molecular gas in ellipticals is not correlated with the old 
stellar population and, hence, is accreted from an external 
source. A similar result for the atomic hydrogen gas was 
reached by Knapp et al. (1985). 

The discovery of kinematically decoupled cores in el- 
liptical galaxies has strongly favoured the idea that these 
systems can be formed through the merging of two disk 
galaxies of similar sizes (Franx & lUingworth 1988, Rix 
& White 1992). The main stellar component is hot and 
slowly rotating while a small centrally concentrated and 
kinematically distinct core exhibits rapid rotation. Numer- 
ical simulations of mergers between two disk galaxies sue- 



cessfully reproduce the observed properties (Hernquist & 
Barnes 1991). The dissipative gas is concentrated in the 
center of the resuhing galaxy and the kinematically de- 
coupled core is observable in both ionized and molecular 
gas (cf. Bcrtola & Bettoni 1988, Balcells & Stanford 1990, 
Wang et al. 1992). The early stages of mergers are ac- 
companied by intense and efficient star formation, giving 
a chaotic appearance to the system. The underlying old 
stellar population, however, has a luminosity profile which 
in most cases closely follows the de Vaucouleurs r^/** law. 
Hence, once the gas has been consumed and star forma- 
tion ceases, the resulting system will not only resemble 
but actually be an elliptical galaxy. 

One elliptical with relatively strong CO emission is 
NGC759 (Wikhnd et al. 1995). We have chosen this 
galaxy for high resolution observations with the aim of re- 
solving the distribution and kinematical structure of the 
molecular gas. In this paper we present our results and dis- 
cuss the rather peculiar conditions under which the molec- 
ular gas exists in this galaxy. 

2. Abell 262 and NGC 759 

NGC 759, a giant EO galaxy, is one of the brightest mem- 
bers of the cluster A262 and forms a noninteracting pair 
with NGC 753. The cluster A262 is comparable to the 
Virgo cluster in its richness and has a similar X-ray lu- 
minosity. The X-ray emission is centered on the elliptical 
NGC 708, which appears to be at rest within the cluster 
(Jones & Forman 1984). NGC 708 is a radio source with a 
dust lane oriented almost perpendicular to its radio axis 
(Ebneter & Balick 1985). Its far-infrared (FIR) luminos- 
ity derived from IRAS data is 3 x 10^ Lq. A deep search 
for CO emission in NGC 708 gave an upper limit to the 
molecular gas mass of 4x 10^ Mq (Braine & Dupraz 1994). 
An optical spectroscopic survey of ten early-type 
galaxies in A262 (including NGC 708 and NGC 759) 
showed that only NGC 759 has a spectrum showing re- 
cent star formation activity (Vigroux et al. 1989). A recent 
study of NGC 759 (E. Davoust private communication) 
shows emission lines of Ha, Nil and SII. The FIR lumi- 
nosity of NGC 759 is 1.1 x 10^° Lq and we observe a CO 
luminosity corresponding to an H2 mass of 2.4 x 10^ M0 
(Wiklind et al. 1995; Tablet). Despite the presence of a 
massive cold ISM component, NGC 759 is a true giant el- 
liptical: It has an r^^^ light profile (i.e. Sandage & Perel- 
muter 1991) and it falls upon the fundamental plane re- 
lationship for elliptical galaxies (cf. Djorgovski & Davis 
1987; Guzman et al. 1993). 

3. Observations and data reduction 

S.l. Observations 

Synthesis observations of NGC 759 were made with the 
Institute de Radio Astronomic Millimctrique (IRAM) in- 
terferometer at Plateau de Bure (PdB), situated at an 



Table 1. Data on NGC 759 



Typo 


EO (El) 




(1) 


Size 


l'.6xl'.4 




(1) 


R.A. (2000) 


0l'^57™50.3= 




(2) 


DEC (2000) 


-|-36°20'35" 




(2) 


Vo 


4665 


kms"^ 


(3.4) 


Distance'^' 


66 


Mpc 


(3) 


R25 


14.5 


kpc 


(4) 


iir-p 


13.24 




(1) 


Mb 


-20.86 




(3) 


{U-B)0, 


1.00 




(1) 


0.52 




(1) 


S12 


< 102 


mjy 


(5) 


S25 


< 132 


mjy 


(5) 


Seo 


732 ± 66 


mJy 


(5) 


Sioo 


2076 ± 187 


mJy 


(5) 


SlAGH^ 


16.4 


mJy 


(6) 


Lb 


3.46 X 10^" 


Lq 


(4) 


ipiR 


1.13 X 10"' 


Lq 


(4) 


^CO 


5.03 X 10* 


Kkms^^kpc"^ 


(4) 


MH2 


2.4 X 10^ 


Mo 


(4) 


Mhi 


< 7.4 X 10* 


Mq 


(7) 


Md 


3.5 X lO*" 


Mq 


(4) 


T/) 


31 


K 


(4) 


SFR") 


7 


Moyr-l 


(3) 


-^pir/Mh2 


4.7 


Lo/Mq 


(4) 


^fir/^b 


0.3 


Lo/iBOL 


(4) 



(1) RC3 (de Vaucouleurs et al. 1991). 

(2) Peak of the CO and radio continuum emission (cf. Fig. 2). 

(3) Adopted in this paper. 

(4) Wikhnd et al. (1995). 

(5) IRAS Faint Source Catalog. 

(6) Feretti & Giovannini (1994). 

(7) Huchtmeier et al. (1995). 

a) Derived from the mean heliocentric velocity of the A262 cluster 
4845 ± 438 (Moss & Dickens 1977) and a Hubble constant of 
75kms~ Mpc~ , corrected for the solar motion relative to the Local 
Group and the peculiar motion of the Local Group with respect to the 
Virgo cluster (Aaronson et al. 1982). 

b) Derived using a v emissivity law. The dust temperature is slightly 
lower than in ref (4) due to modified FIR fluxes. 

c) Derived using SFR ^ 6.5 X 10"^" LpiR Mq yr"^ 



altitude of 2560 m. The PdB array comprises 4 telescopes 
of 15-m diameter which can be placed at 26 different sta- 
tions along a T-shaped runway extending 288 m East- 
West and 160 m North-South (for a detailed description, 
see Guilloteau et al. 1992). NGC 759 was observed with 
the standard CD configuration (see Table 2). Maximum 
and minimum baselines were 160 and 24 m, respectively. 
The observations were made on 7 occasions between June 
1993 and March 1994, with a total observing time of 52.5 
hours. The SIS receivers were single-sideband tuned to the 
redshifted ^^CO, J=l-0 frequency 113.508 GHz, with the 
CO line placed in the upper sideband. The receiver tem- 
peratures were 40-140 K, resulting in system temperatures 
of 300-800 K (T^). The correlator was configured as five 
slightly overlapping 64 x 2.5 MHz subbands, giving a ve- 



36 20 40 



36 20 20 



264.1 Q - ^ 




1 57 51 57 5Q 57 49 



-200 200 



Fig. 1. Channel maps of the CO emission in NGC 759, separated by lOMHz (26.4knis~^). Velocity offsets are relative to the 
systemic velocity of 4665 km s~^. The cross in each map corresponds to the phase reference center. The contour levels are in 
steps of lOmJy (0.13 K). Also shown are the dirty beam and the uv-coverage. The latter is in units of meters. Coordinates are 
J2000.0. 



locity resolution of 6.6 km s and a total bandwidth of 
1340kms"\ 



3.2. Data reduction 

The CLIC programme and associated Grenoble software 
were used for all stages of the data reduction. Bad 
scans were flagged, followed by bandpass calibration us- 
ing strong continuum sources observed in the beginning 
of each observing run. Each correlator subband was fitted 
separately with a 7-lOth order polynomial. The edges and 
two channels in the center of the subband were affected 
by the Gibbs phenomenon and were omitted in the fitting 
process. The errors in the amplitude and phases for the 



bandpass correlations were 2 — 8% and 1.2 — 4°, respec- 
tively. 

Cubic splines were fitted to the measured phases of the 
phase calibrators. In a few cases a focus measurement dur- 
ing the observations introduced a phase jump. In order to 
match these discontinuities we applied two fits separately. 
Due to bad weather conditions the observations in June 
27, August 31, and September 3 had residual phase er- 
rors of 25 — 35°. The rest of the observations have residual 
phase errors of 5 — 15°. 

Long integrations on cither 3C345, 3C454.3 or 3C84 
were made at the beginning of each observing run. These 
sources were used as bandpass calibrators as well as flux 
references. Since most calibrator sources at millimeter 
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Fig. 2. Total integrated CO(l-O) emission obtained by adding all channels between — 237knis^^ and +185 kms^^ (relative to 
the systemic velocity). The level spacing is O.SJykms^^ and the rms noise in the map is 0.8 Jykms^^. The cross marks the 
phase reference center and the optical position, but is offset a few arcsec from the CO and radio continuum peaks (see text for 
details). Coordinates are J2000.0. 



wavelengths are variable, it is difficult to apply a correct 
flux scale. The flux values were bootstrapped from mea- 
surements of sources with accurately known fluxes. For 
3C84 we used 5.1±0.8 Jy and for 3C454.3 6.9±0.8 Jy, de- 
rived with the IRAM 30-m telescope in August 1993. We 
could not get a reasonable solution using 6.2 Jy for 3C345, 
derived in February 1993. Instead we used 0235-1-285, with 
2.05 Jy derived in July 1993. For the March 29 1994 obser- 
vation 0133-1-476 was used as a flux calibrator (1.45 Jy). 
The adopted fluxes were used to calculate the antenna 
gains. With these gains the fluxes of all phase calibrators 
were determined, only using scans with the highest ampli- 
tudes to ensure that no phase fluctuations had reduced the 



amplitude. Cubic splines were then fltted to the NGC 759 
data. 

From the calibrated data we constructed uv-tah\es, ap- 
plied a fast Fourier transformation and deconvolved the 
resulting maps with the CLEAN algorithm using nat- 
ural weighting. The resulting size of the clean beam is 
3'.'lx2'.'3, corresponding to 990 x 735 pc at a distance of 
66 Mpc. The dirty beam and the wt^coverage are shown 
as part of Fig. 1. 



4. Results 



Table 2. Observing parameters 



Date 


Configuration"' 


Hours 


Bandpass cafibrator 


Phase cafibrator 


27 June 1993 


N05 N13 ElO 


8.5 


3C345 


0234+285 


1-2 Jul 1993 


N05 N13 ElO W12 


9.2 


3C454.3 


0234+285 


19-20 Aug 1993 


E03 W05 N05 


7.4 


3C84 


3C84 


21 Aug 1993 


E03 N05 WOO 


5.0 


3C84 


3C84 


31 Aug - 1 Sep 1993 


ElO W09 N15 N03 


10.7 


3C84 


3C84 


3-4 Sep 1993 


ElO W09 N15 N03 


9.6 


3C84 


3C84 


29 Mar 1994 


E03 N05 WOO W05 


2.2 


3C345 


0133+476, 3C84 



a) N, W and E refer to the northern, western and eastern arm, respectively. Numbers denote 
the distance from the intersection of the arms (WOO) in units of 8 m. 



4-.1. The CO distribution 

Individual channel maps with a velocity separation of 
26.4kms~^ (lOMHz) are shown in Fig. 1. CO emission 
is clearly detected in channels approximately ±200 km s^ 
around the systemic velocity of 4665 km s~ . The strongest 
emission is seen around +185 and — 185kms~ relative to 
the systemic velocity, with peak flux densities of 95 and 
65 mjy per beam. The rms noise in the individual maps 
is 6 mJy per beam. 

A map of the total integrated CO emission was 
made by adding all the channels between — 237kms~ and 
+ 185kms~ . The result is presented in Fig. 2. The peak 
of the CO emission is 9.2 Jykms"^. By summing the same 
number of channels without emission we obtained an rms 
noise of 0.8 Jykms"^. 

Two things are immediately evident from Fig. 2. First 
of all, the peak of the CO emission is offset ^ 4" from the 
nominal pointing of the telescope. Since the peak of the 
CO emission also corresponds to the kinematical center, as 
well as the center of radio continuum emission (Feretti & 
Giovannini 1994), this offset most likely reflects an uncer- 
tainty in the optical coordinates. Secondly, the emission 
region has an elliptical shape, with an extent of 5 '.'4x4 '.'3 
(1.8 X 1.4 kpc) at the 50% level. This is larger than the 
synthesized beam which has a FWHM size of 3'.'lx2'.'3. 
However, the synthesized beam and the CO emission re- 
gion have similar ratios of the major and minor axis as 
well as the same position angle (~ 50°). This could be a 
coincidence or an indication that our synthesized beam is 
larger than estimated. In addition to the bright compact 
component, there is low level emission extending towards 
the South- West, to a distance of ^ 7" from the center. 

At a frequency of 113.5 GHz and with the CD config- 
uration, the PdB interferometer will not detect emission 
more extended than about 30", corresponding to ^ 10 kpc 
at a distance of 66 Mpc. We can check if such extended 
emission exists in NGC 759 by comparing our single dish 



flux with the flux obtained with the interferometer. Inte- 
grating the emission above 0.8 Jykms^^ for the interfer- 
ometer data, we obtain a total flux of 37±2 Jykms^^. The 
single dish CO (1-0) spectrum obtained with the IRAM 
30-m telescope has an integrated flux of 11.8Kkms~^ 
(Wiklind et al. 1995). Converting the single dish CO(l-O) 
flux to Jy using 4.47 Jy/K, we get 53 ± 3 Jy. The inter- 
ferometer thus recovers ^ 70% of the single dish flux. It 
is possible that 30% of the emission seen with the sin- 
gle dish telescope is uniformly extended over more than 
10 kpc; however, given the uncertainties associated with 
calibration (about 30% for both single dish and interfer- 
ometric data) and the fact that the emission profiles are 
very similar in shape (see Fig. 3), we believe that it is more 
likely that almost all molecular gas seen with the single 
dish telescope is concentrated within the central 5". 

4-2. The velocity field 

The channel maps in Fig. 1 show strong CO emission in 
channels at approximately ±185 km s^ from the center 
velocity. These channels correspond to the two peaks in 
the double-horned profile shown in Fig. 3a. A close in- 
spection of the channel maps shows that the two main 
emission regions are separated by a few arcseconds. The 
separation is, however, of the order of the angular reso- 
lution. We therefore made Gaussian fits to the individual 
channels, finding a systematic shift in center velocity with 
position, with redshifted emission coming from the North- 
East and blueshifted from the South- West (Fig. 4) . From 
a fit to the velocities in Fig. 4 we derive a position angle 
of ^ 50°, measured from North to East. 

A position-velocity diagram along a position angle of 
50° is shown in Fig. 5, where the offsets in position are 
relative to the peak of the integrated emission. The two 
peaks of CO emission seen in Fig. 5 are symmetrical with 
respect to the center, both in position and velocity. The 
separation of the two peaks is 4", corresponding to 1.3kpc. 
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Fig. 3. a) CO(l-O) spectrum integrated over the emission re- 
gion seen in Fig. 2. b) Single dish CO(l-O) spectrum obtained 
with with the IRAM 30-m telescope, c) Single dish CO(2-l) 
spectrum obtained with the IRAM 30-m telescope. 



Although the extent is uncertain due to possible residual 
phase errors, which tend to enlarge the emission region, 
we will use 2'.'0 (650 pc) as an 'effective radius' of the 
molecular gas distribution. This agrees with the average 
radius of a model fitted to the data (see Sect. 5.3). 

4-. 3. The CO J=2-l/J=l-0 line intensity ratio 

The extent of the CO emission region of about 4-5" en- 
ables us to derive a size-corrected CO J —2—1 / J —\—Q line 
intensity ratio. The single dish CO J=l— and J=2— 1 
spectra have integrated line intensities of 11.8 K km s^^ 
and 15.2Kkms~ , respectively. Convolving the Gaussian 
telescope HPBW of 23" for the CO(l-O) line and 12" for 
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Fig. 4. Positions of Gaussian fits to the individual 10 MHz 
channel maps relative to the phase reference center. The trian- 
gle represents the peak of the integrated CO emission seen in 
Fig. 2 and coincides with the kinematical center. The line is a 
best fit position angle to the kinematical data. 
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Fig. 5. Position-velocity diagram along the major axis indi- 
cated in Fig. 4. The level spacing is 8mJy/beam. The offsets 
are measured relative to the position of maximum integrated 
emission (cf. Fig. 2). Velocities are relative to the systemic ve- 
locity. 



the CO (2-1) line with a Gaussian shaped source distribu- 
tion of HPBW 4"0, we get a J=2-l/J=l-0 line ratio 
of 0.4. This value is much lower than what is expected 
from optically thick and thcrmalized CO emission and, 
since the source size is small compared to both telescope 
beams, represents a good estimate. 
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Fig. 6. The rotation curve corresponding to a spherical galaxy 
with an r^^'* projected luminosity profile (Sect. 5.1). The an- 
gular scale on top corresponds to a distance of 66Mpc. 
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4-.4- Comparison with radio continuum 

The 1.4 GHz radio continuum emission of NGC 759 has 
been resolved with the Very Large Array (Feretti & Gio- 
vannini 1994). Its structure is characterized by diffuse 
emission without a compact core. The ratio between radio 
and FIR fluxes is in agreement with that of star forming 
galaxies. A comparison of the distribution of the radio 
continuum emission and the molecular gas shows that the 
continuum emission is centered on the peak of the molecu- 
lar gas distribution and that it is contained within the CO 
emission region (see Fig. 1 in Feretti & Giovannini 1994). 
If the continuum emission at 1.4 GHz were due to 
free-free emission and the frequency turn-over lower than 
1.4 GHz, we would expect a flux density of ^ lOmJy 
at 113.5 GHz due to thermal emission in HII regions. 
A higher turn-over frequency would yield a higher flux 
density. We did not detect any continuum emission at 
113.5 GHz with a 3a limit of 3.9mJy/beam. It is there- 
fore likely that the continuum emission at 1.4 GHz is dom- 
inated by synchrotron radiation. 

5. Dynamics and kinematics of NGC 759 

5.1. A mass model 

In order to interpret our CO data, we need to construct 
a reasonable model of the mass distribution in NGC 759. 
Since it is classified as an EO galaxy, a spherical model with 
a light distribution following the de Vaucouleurs r^'^ law 
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Fig. 7. Synthetic spectra of an unresolved source, derived using 
the formalism presented in Sect. 5.3. a) A filled disk and, b) a 
ring. Only the ring is compatible with the observed emission 
profile from NGC 759. The numbers in the spectra correspond 
to the inclination of the gas distribution. 



is a reasonable choice. We will use the potential-density 
pair presented by Hernquist (1990), which closely approx- 
imates the de Vaucouleurs r^^^ law for elliptical galax- 
ies. This particular model has the advantage of expressing 
most aspects of the dynamics in an analytical form. 
The mass as a function of radius is given by 



M{r) 



Mr^ 



(r + a)2 



(1) 



where a is a parameter related to the effective radius 
To within which half of the luminosity is confined, r^. w 
1.8153a. M is the total mass, which we derive by first cal- 
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Fig. 8. Synthetic spectra obtained from the model gas distri- 
bution of NGC 759 (Sect. 5.3), viewed at different inchnations. 



culating the global M/L ratio. Following Binney (1982), 
we have for a spherical r^'^ galaxy, 



M 



0.047ct2 

I a Re. 



(2) 



where i?e is the effective radius in kpc , Ig is the surface 
brightness at Re in units of Lq pc~^ and a is the stel- 
lar velocity dispersion in knis^ . The corresponding val- 
ues for NGC 759 are R^ = 2.6 kpc and /„ = 225L0pc-2 
(Sandage & Perelmuter 1991). The velocity dispersion was 
estimated by Vigroux et al. (1989) to be ct = 250 kms"^. 
Since their value could be influenced by a higher than av- 
erage velocity dispersion in the center of NGC 759, we will 
use a = 225 km s^^ as a mean value. This gives M/Lb ~ 4 
in solar units and a mass of 1.4 x 10"'^"'^ Mq within a radius 
corresponding to R25 (14.5 kpc). Extrapolating to i?oo, the 
total mass becomes 1.7 x 10^^ Mq. The circular velocity 
as a function of radius is 



Kot, - 



VGMr 



The rotation curve for NGC 759, using Eq. (3), is shown 
in Fig. 6. The maximum rotational velocity of 360 km s^ 
is reached at a galactocentric radius of 1.4 kpc. The rota- 
tional velocity at r = 650 pc, the assumed extent of the 
molecular gas, is 330 km s" . Assuming that the molecular 
gas is distributed in a disk or ring -structure (see Sect. 5.3) 
we derive an inclination of 40° from the width of the ob- 
served emission profile. 



5.2. The gas mass fraction 

The total mass enclosed within a radius of 650 pc (Mdyn) 
is given by Eq. (1) as 1.7 x 10^" M©. This constitutes 10% 
of the total galactic mass and results from the assumption 
of a constant M/L ratio and the strongly peaked luminos- 
ity distribution of elliptical galaxies. For elliptical galax- 
ies in general, the M/L ratio tends indeed to be constant 
within ~ 1.2 i?c and near a value of 4Mo/Lq (Bertola et 
al. 1993). The Mn^/Mdyn ratio becomes 18% within the 
central 650 pc, where we have included helium at a mass 
fraction of 1.3. As long as the mass model for NGC 759 
gives the correct dynamical mass, the ratio Mn^/M^yn 
only depends on the radius of the molecular gas distri- 
bution. Whereas the gas mass fraction in the center of 
NGC 759 is comparable within a factor two to the gas 
mass fractions found in normal spiral galaxies, it falls to 
less than 2% when including the entire galaxy. The mass 
of atomic hydrogen is < 7.4 x 10^ Mq (Huchtmeier et al. 
1995) and cannot greatly influence the Mgas/Mdyn ratio. 
Hence, NGC 759 is a gas-poor galaxy when viewed on a 
global scale. 

The H2 mass of NGC 759 is derived using a standard 
conversion factor between integrated CO emission and the 
column density of H2 (cf. Wiklind et al. 1995). This ratio 
is unlikely to be universal, depending on excitation condi- 
tions and metallicity. For instance, if the metallicity of the 
molecular gas in NGC 759 is low compared to our Galaxy, 
we may underestimate the mass by a factor of ^4. This 
uncertainty in the molecular gas mass does not greatly 
influence the results discussed in this paper. The same 
uncertainty also applies to other galaxies with which we 
compare NGC 759. 

5.S. The CO line profile and the molecular gas distribution 

The double-horned CO profile of NGC 759 is already in 
itself a clue to the molecular gas distribution, ruling out 
too large gas concentrations in the very center. It can be 
shown by simple integration of the emission along the line 
of sight that the bulk of the CO gas should lay outside 
the rigid part of the rotation curve to produce the double- 
horned profile. For an assumed axisymmetric gas distri- 
bution of surface density n{r), a typical global spectrum 
dN/dv{V) is the sum over all radii r of 



(3) dN 



dV 



n(r)rdrdd = / n{r)rdr 



dV 



(4) 



dv J J V^ot {r) sin 9 sin i 

where V,:ot{'>') is the rotational velocity at radius r, and i 
the inclination of the galaxy in the sky (i = is face-on) . 
Since the observed V = K-ot cos 9 sin i, it is straightforward 
to derive that in the simple but unphysical case of a con- 
stant rotational velocity VJ-ot at all radii, the double horn 
profile 

-1/2 

1 - ( T^-^ 1 , (5) 



dN 



dv 



K-ot sin i 



results whatever the radial distribution of the surface den- 
sity n{r). In the opposite case of a rigid body rotation 
curve, and with a constant gas surface density n{r) — uq 



until Rn 



where the rotation curve turns over to become 



flat, the spectrum becomes 



dv 



1- 



V 



K. 



sini 



-,1/2 



(6) 



where l^nax is the rotational velocity reached at i?max- 
These cases represents two extremes of rotation curves and 
show that a parabolic profile (large i?max) will transform 
into a double-horned when i?niax decreases. The depth of 
the double-horn profile is thus a good indicator of the gas 
concentration with respect to the turn-over of the rotation 
curve. Emission profiles obtained from Eqns. (5) and (6) 
are shown in Fig. 7. 

We have modeled in more detail the gas distribu- 
tion and kinematics to better determine the constraints 
brought by the CO observations. We have assumed that 
the gas is in cylindrical rotation in a thick disk or ring, in 
the potential of the elliptical galaxy, as fixed in Sect. 5.1, 
also accounting for the potential of the gas itself. The lat- 
ter is not a negligible contribution, given its high cen- 
tral condensation. For the disk we choose an analytical 
density-potential pair in the form of Toomre disks of or- 
der n=2 (cf Toomre 1964). To simulate a disk, we choose 
the surface density: 



n(r) = "0 ( 1 + "H 



-5/2 



(7) 



where d is the disk characteristic scale, and the corre- 
sponding rotational velocity is 



V'{r) = 



SttGuq 
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(8) 



For a ring, we choose the difference of two such distri- 
butions with di and d2 as characteristic scales, and with 
the same central surface density uq. Since the total gas 
mass is fixed (Mh2 = 2.4 x 10^ Mq), and so is the ellip- 
tical potential, we only vary the scales di and ^2- The 
global spectrum is then derived by a simple radial inte- 
gration. A best fit with the observed spectrum is obtained 
for a ring constructed with scale lenghts di = 800 pc and 
^2 = 750 pc, viewed at an inclination i — 40° (see Fig. 8). 
This gives a maximum surface density of 750 Mq pc~^ at 
r — 500 pc. The surface density profile is shown in Fig. 9. 
Although the ring has a maximum at 500 pc, the distribu- 
tion is extended at galactocentric distances up to ^1.5 kpc 
and the radius of 650 pc inferred from the interferometer 
data represents an average radius for the gas distribution. 
Is this extended gas distribution in agreement with the 
barely resolved interferometric map? We have answered 
the question by simulating the ring with TV =10 000 par- 
ticles, distributed as above, with a cylindrical rotation in a 



disk thickened by a sech^(z/zo) law, with zq — 200 pc, and 
with an isotropic velocity dispersion of 30 km s^ . We then 
'observed' the distribution with an inclination i = 40°, 
and an elliptical beam of 3'.'lx2'.'3. The total integrated 
density of the model distribution is shown in Fig. 10, and 
is in good agreement with the observations. 




Radius (kpc) 

Fig. 9. The gas surface density as a function of galactocentric 
radius, obtained using the ring constructed from two Toomre 
disks (cf. Sect. 5.3). 
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Fig. 10. A synthetic image of the ring-structure seen at an 
inclination of 40° , convolved with a Gaussian beam of the same 
HPBW size as the CLEANed beam (S'.'l x 2 "3). 



Table 3. H2 densities and rB(l-O) for difTerent fractional CO abundances and for CO J=2— 1/J=1— = 0.4 
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a) Tb is the excess brightness temperature above the microwave background. 

b) M^^/Lqo = 2.1^Jn{n2)/T^{l - 0) (ef. Radford et al. 1991). 
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Fig. 11. The ratio of the molecular gas surface density Sgas 
(including helium at a mass fraction of 1.3) and the critical gas 
surface density Ecrit, as a function of galactocentric radius for 
a galaxy described by the mass model presented in Sect. 5.1. 



6. Discussion 

6.1. Is NGC759 a starburst elliptical? 

Apart from the mere presence of molecular gas, there 
are three facts which suggest that massive star forma- 
tion is taking place in the inner 650 pc of NGC759: (i) 
optical emission lines from HII regions, (ii) diffuse radio 
continuum emission and (Hi) FIR emission which most 
likely comes from dust heated by massive stars. With a 



Salpeter initial mass function, the inferred star formation 
rate (SFR) is ~7MQyr-i (Tabid). 

How does the molecular data fit into the picture of 
NGC 759 as an actively star forming elliptical galaxy? Be- 
low we discuss the molecular properties in more detail in 
view of its implications for the star formation activity in 
NGC 759. 

6.2. The molecular gas surface density and Ecrit 

Kennicutt (1989) has argued for the existence of massive 
star formation activity in spiral galaxies only when the 
gas surface density exceeds a critical value, for which the 
gaseous disk is expected to become unstable to large scale 
density perturbations. Using the Toomre (1964) formula- 
tion, the critical surface density can be expressed as: 



ttG" 



(9) 



where a is a constant of order unity (0.7 for spiral galax- 
ies), K is the epicyclic frequency, c is the gas velocity dis- 
persion and G the gravitational constant. Using the rota- 
tion curve corresponding to the mass model presented in 
Sect. 5.1, the critical gas surface density becomes 
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.Vrr 
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0.55rc 
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(10) 



where r is the galactocentric distance expressed in kpc, 
VJ-ot is the rotational velocity in kms^ , and we have as- 
sumed a gas velocity dispersion of 30kms^ . The ratio of 
the molecular gas surface density obtained from the model 
distribution presented in Sect. 5.3 and the Ecrit derived us- 
ing Eq. (10) is shown in Fig. 11. Due to the rapid increase 
of Ecrit with decreasing radius, the ratio is everywhere less 
than one. Despite the high molecular gas surface density 
we do not expect the gas to be highly susceptible to large 
scale density perturbations. The ratio is, however, of or- 
der unity around the peak of the molecular ring, and this 



allows gravitational instabilities to develop on the same 
scale as in the disk of normal spiral galaxies. The maxi- 
mum of Sgas/Scrit fs-Us at a slightly larger galactocentric 
distance than the peak of the molecular ring. A thin ring 
of HII regions should therefore be visible in narrow band 
Ha images at a distance of ^2'.'0 from the nucleus of 
NGC759. 

6.3. Eh2 o,^d ih^ L-pm/MYi^ ratio 

The molecular gas surface density of 750 M© pc^^ found 
for NGC 759 is high when compared with the Galactic 
molecular ring (12M0pc^^) and the Galactic center re- 
gion (200 - 300 M© pc"2), but it is a factor 2-3 smaller 
than the surface density obtained by assuming a constant 
surface density disk truncated at a radius 650 pc. In the 
latter case NGC 759 would have a gas surface density sim- 
ilar to those of ultraluminous FIR galaxies (cf . Scoville et 
al. 1991). The difference comes from the larger extension of 
the CO distribution at low surface density (out to 1.5 kpc) 
in the modelled ring distribution (see Sect. 5.3). 

In Fig. 12 we compare the Lfir/Mh2 ratio vs. Sua for 
ultraluminous FIR galaxies and NGC 759, with Eh2 de- 
rived assuming a truncated gas disk and using the virial 
theorem. NGC 759 has an Lfir/Mh^ ratio of 4.7 Lq/Mq, 
which is much smaller than for ultraluminous FIR galax- 
ies who typically have ratios in excess of 40 Lq/Mq. In 
fact, the Lfir/Mh2 ratio of NGC 759 is only half the av- 
erage value for normal IRAS selected spiral galaxies (cf. 
Young et al. 1989). The large Sh2 of NGC 759 and the 
low Lfir/Mh2 ratio suggests that NGC 759 is undcrlu- 
minous in Lfir with respect to its H2 mass as compared 
to ultraluminous FIR galaxies. This could be caused by a 
lack of massive stars heating the dust in NGC 759, sug- 
gesting a low star formation efhciency, or the presence of 
an extra heating source of the dust in ultraluminous FIR 
galaxies, possibly an AGN. If the metallicity of the molec- 
ular gas in NGC 759 should be very low, which it could 
be if the gas was accreted from a dwarf galaxy, the H2 
mass would be underestimated, enhancing the difference 
between NGC 759 and the ultraluminous FIR galaxies. 

6.4^. Subthermal excitation of the CO emission 

A CO J=2-l/J=l-0 line ratio as low as 0.4 (Sect. 4.3) 
is rarely encountered in molecular clouds, neither in our 
Galaxy nor in external spiral galaxies (Braine & Combes 
1992). The low ratio implies a subthermal excitation of 
the J= 2 level and is usually interpreted as cold and/or 
diffuse molecular gas. 

Table 3 summarizes large velocity gradient calculations 
for a CO J=2-l/J=l-0 ratio of 0.4. The calculations 
are taken from Castets et al. (1990) and give for three 
different fractional CO abundances per velocity interval, 
X (CO) /{dV/dr) (kms^^ pc"^)""'^, the brightness temper- 
ature Tb, the II2 density and the ratio Ma^/Lco- -'^(CO) 



is equal to n(C0)/n(Il2). A small fractional CO abun- 
dance can be interpreted as a low metallicity gas and a 
large abundance as a high metallicity gas. The M^i^/Lco 
ratio is derived in the same way as in Radford et al. 
(1991). A value of 4.5Mq/Lq corresponds to the stan- 
dard Nii.^/Ico conversion ratio. 

A fractional CO abundance lower 

than 10^^ (kms^^ pc^^)^^ can be excluded on grounds 
of a very low brightness temperature and an Mh^^/Lcq 
ratio ^10"^ Mq/Lq. In the case of a small fractional CO 
abundance we underestimate the molecular mass by a fac- 
tor 5-10, whereas a high fractional abundance gives the 
correct mass to within a factor of 2. In the simple hy- 
pothesis of a homogeneous one-component medium, un- 
der the LVG approximation, two extreme physical con- 
ditions are possible for the molecular gas: either cold and 
dense (Tk « 10 K, n(H2) <; 10'^ cm""^) or warm and diffuse 
(Tk w 40K, n(Il2) w 10^ cm^-^). A more realistic solution 
might be a multi-medium composed of dense clumps em- 
bedded in a diffuse component, where the CO emission 
is dominated by the diffuse component while most of the 
mass is contained in the dense clumps. This can be seen 
from the higher M-h^/Lfi-r, ratio of the cold and dense 
component compared to the diffuse gas (see Table 3). 

6.5. Comparing with a bonafide merger: NCC 7252 

A good example of an advanced merger is NGC 7252, 
which shows evidence of a recent starburst in its center (cf. 
Schweizer 1982). Although the merger initially consisted 
of two disk galaxies, the overall light distribution now fol- 
lows the r^/^ law (Schweizer 1982, Stanford & Bushouse 
1991). 

Dupraz et al. (1990) observed NGC 7252 in the CO 
J=l— and J=2— 1 lines with a single dish telescope 
(SEST), and Wang et al. (1992) observed it with the 
Owens Valley interferometer, achieving a resolution of 
3'.'8x2'.'l (1.1 X 0.6kpc, assuming a distance of 63Mpc). 
The molecular gas distribution in NGC 7252 has a strik- 
ing resemblance to that of NGC 759. The CO profiles from 
the single dish and the interferometer are double horned, 
and the inclination of the gas appears to be around 40° . 
At least 72% of the CO(1~0) emission in NGC 7252 comes 
from within a radius of 1.5 kpc. The CO J=2— 1/J=1— 
ratio, corrected for the size of the emission region, is 
0.4 — 0.6 (depending on how much of the C0(2~l) emis- 
sion that is extended). The molecular gas surface density 
of NGC 7252 is similar to NGC 759, when calculated in 
the simple way discussed in Sect. 5.2. However, while the 
dust temperatures are comparable, the Lpi^/ Mn^ ratio of 
NGC 7252 is a factor 3 higher than in NGC 759 

Since the projected light distribution of NGC 7252 
closely follows the r^^^ law, it is appropriate to use a mass 
model similar to the one we used for NGC 759 (Sect. 5.1). 
From the light profile presented by Stanford & Bushouse 
(1991) we derive an effective radius Ve of lO'.'S (3.3 kpc). 



If we express the total mass of NGC 7252 as a function of 
the mass-to-hght ratio, 



Mt. 



^h' 



(11) 



and use Eq. (1) and (2), wc get the following expressions 
for the mass within radius r = 1.5 kpc and the correspond- 
ing rotational velocity, 



M{r) = 9.4x 10'''( — ) M0, 



V{r) = 166^ 






kms 



(12) 



(13) 



The H2 mass within a radius of 1.5 kpc is 3 x 10® M0 
(Wang et al. 1992), giving a A^HaZ-^dyn ratio of 
0.3 (M/Lb)~ • The mass-to-light ratio must be at least 3 
to give the observed velocity width {i = 90°). This gives a 
molecular gas mass fraction of ^ 10% in the inner region 
of NGC 7252. This is considerably smaller than the 45% 
derived by Wang et al. (1992). If the inclination is in the 
range 40° — 50°, as implied by the interferometer map, 
the mass-to-light ratio is 5 — 7, further diminishing the 
molecular gas mass fraction. 

Apart from the Lfir/-'^H2 ratio, the advanced merger 
NGC 7252 and the elhptical NGC 759 show striking simi- 
larities in their molecular gas properties. 

The very high gas mass fractions found in the central 
regions of ultraluminous FIR galaxies (cf. Scoville et al. 
1991), including NGC 7252, have been obtained from a 
dynamical mass derived using the virial theorem. This is 
applicable to disk systems with a small central mass con- 
centration compared to elliptical galaxies. However, most 
of these ultraluminous FIR galaxies are mergers, and pre- 
sumably have a light distribution following the r^'"* law. 
Hence, it would be appropriate to use a more concentrated 
mass model such as that described by Hernquist (1990). 
This would decrease the central gas mass fractions in ul- 
traluminous FIR galaxies by a factor 3 — 10. 

7. Summary and conclusions 

NGC 759 is an elliptical galaxy which contains 2.4 x 
10® Mq of molecular gas, most of it confined to a cir- 
cumnuclear ring with an average radius of ~^650pc. The 
average molecular gas surface density in the center re- 
gion of NGC 759 is 750Mopc~2. Although this value is 
high, it is always less than the critical gas surface den- 
sity for large scale gravitational instabilities. The low 
CO J=2— 1/J=1— line ratio of 0.4, is consistent with a 
molecular gas consisting of cold and dense clumps (Tk ~ 
lOK, n(H2) ^ 10'^ cm^"^) embedded in a warm and dif- 
fuse molecular gas (Tk ~ 40 K, n(H2) « 10^ cm"^). The 
dense component contains most of the molecular gas mass. 
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Fig. 12. The Lfir/Mhj ratio as a function of the aver- 
age nuclear molecular gas surface density for ultraluminous 
far-infrared galaxies (open circles), NGC 7252 (filled circle) 
and NGC 759 (star). The gas surface density has been derived 
in the same manner for all galaxies; using a truncated gas disk 
and the virial theorem. 



while the diffuse gas dominates the emission. Finally, when 
compared with galaxies of similar gas surface densities, 
NGC 759 appears to be underluminous in ipiR with re- 
spect to its molecular gas mass by a factor 3 — 10. 

In a quiet environment, molecular gas condenses into 
a very clumpy structure with local high densities and a 
low volume filling factor, resulting in a low average gas 
surface density. Intense star formation partially destroys 
the clumpy structure and the fraction of diffuse molecular 
gas is considerably enhanced. Consequently the filling fac- 
tor and the average surface density will also increase. Al- 
though the kinetic temperature increases due to these pro- 
cesses, the CO excitation temperature can become quite 
low due to subthermal excitation caused by the low vol- 
ume density of H2 . 

This scenario, which seems to be applicable to 
NGC 759, implies that star formation is ongoing and has 
been so long enough to influence the molecular gas prop- 
erties in the circumnuclear ring. The low Lfir/-^H2 ratio 
could have its explanation in the lack of an additional 
heating source for the dust in the form of an AGN. 

The high concentration of molecular gas in the center 
of NGC 759, the lack of atomic gas and the relatively mod- 
est amount of star forming activity is reminiscent of the 
conditions found in the merger NGC 7252. The main dif- 
ference being that while NGC 7252 exhibits gas-rich tidal 



tails and a high star formation efficiency, NGC 759 shows 
no sign of gravitationally induced disturbances and has a 
low star formation efficiency. 

The cluster A262 contains a high fraction of spiral 
galaxies (75%). The spirals are less centrally concentrated 
and have a higher line~of-sight velocity dispersion then 
the eUipticals (Moss & Dickens 1977). NGC 759 is situated 
in outer parts of A262, at a projected distance of 1.4Mpc 
from the cluster center. Could NGC 759 represent a very 
late stage in the merging of two disk galaxies? The fate 
of the gaseous component in the merging of two gas rich 
galaxies is two-fold: a major part of the gas looses angu- 
lar momentum quickly and falls towards the center while 
the merging process is on-going. If this gas component 
is capable of forming stars, it can produce the high stel- 
lar phase-space density characteristic for elliptical galax- 
ies. The other part of the gas is ejected into tidal arms, 
which remain bound to the system (Hernquist & Spergel 
1992; Hibbard & Mihos 1995). This gas will fall back to 
the galaxy at later times. If elliptical galaxies are formed 
through the merging of two gas rich disk galaxies, an im- 
portant question is: where are the galaxies which are in 
between the final product, an inert and gas-free elliptical 
galaxy, and those which still show signs of the merging 
process, such as NGC 7252. This question was first raised 
by I. King and has been named the 'King Gap'. Galaxies 
in the 'King Gap' should look like proverbial elliptical but 
still retain some signatures of the merging process. These 
signatures can be: 

— Gas and dust components as well as an intermediate 
aged stellar population. 

— Optical fine structure in the form of faint shells, ripples 
and boxy isophotes (cf. Schwcizer & Seitzer 1992). 

— Kinematically distinct sub-systems (cf. Surma & Ben- 
der 1995). 

— Gas in the outer part of the tidal arms require very 
long time scales to return to the galaxy. Hence, one 
would expect the presence of low HI column density 
at large galactocentric distances (cf. Hibbard & Mihos 
1995). 

NGC 759 does show unusual gas properties for an el- 
liptical galaxy and if it represents a galaxy in the 'King 
Gap', one or more of the additional signatures should be 
found in a more extended study. 
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